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Abstract. The Brieskorn lattice H" of an isolated hypersurface 
singularity with Milnor number /i is a free (D{{s}}-niodule of rank 
fjL with a differential operator t = s^dg- Based on the mixed Hodge 
structure on the cohomology of the Milnor fibre, M. Saito con- 
structed C{{s}}-bases of H" for which the matrix of t has the form 
A = Ao + Ais . We describe an algorithm to compute the matri- 
ces Aq and Ai. They determine the differential structure of the 
Brieskorn lattice, the spectral pairs and Hodge numbers, and the 
complex monodromy of the singularity. 



1. The Milnor Fibration 

Let / : (C""'"^,0) — > (C,0) be a holomorphic function germ with 
an isolated critical point and Milnor number ^ = dim^ <C{x} / {d{f)) 
where complex coordinate system of (C^^^jO) and 

d. = dxQi ■ ■ ■ 1 dx„- By the finite determinacy theorem, we may assume 
that / G C[s]. By E.J.N Looijenga [71 2.B], for a good representative 
/ : X — y T where T C C is an open disk at the origin, the restriction 
f:X' — toT' = T\{0} and X' = X\f-\0) is a fibre bundle 
unique up to diffeomorphism, the Milnor fibration. By J. Milnor [01 
6.5], the general fibre Xt = f~^(t), t G T', is homotopy equivalent to 
a bouquet of /i n-spheres and, in particular, its reduced cohomology 
is Yl^{Xt) = 5fc,n^^ where S is the Kronecker symbol. Since T' is 
locally contractible, the n-th cohomologies H(t/) = H"(Xt/) of Xu = 
f~^{U) form a locally free Z-sheaf of rank n and He = H 0^ C is a 
complex local system of dimension fi. Hence, the sheaf of holomorphic 
sections = H (g)^ of He is a locally free ^T'-sheaf of rank /i, 
the cohomology bundle. By P. Deligne [H 2.23], there is a natural flat 
connection V : Jif — > M' fi^, on ^ with sheaf of fiat sections 
H = ker(V), the Gauss-Manin connection. 



2. The Monodromy Representation 

Let t be a complex coordinate of T C C, i : T' — > T the canoni- 
cal inclusion, and u : T°° — > T' the universal covering of T' defined 
by u{t) = exp(27rir) for a complex coordinate r of T°° C C. Then 
the covariant derivative Vot of V along dt induces a differential op- 
erator dt on i,,J^ and the pullback f°° : X°° = X' Xt' — > 
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is a fibre bundle with = X^^r)-, the (canonical) Milnor fibre. 
Since T°° is contractible, the n-th cohomologies H{U) = H"(X^) of 
= {f°°)~^{U) form a free Z-sheaf of rank /z and u^:H is the sheaf 
of multivalued sections of H. Lifting closed paths in T' along sections 
of H defines the monodromy representation 7ri(T',t) — > Aut(Ht) on 
Ht inducing the monodromy representation 7ri(T') — ^ Aut(if) on the 
cohomology H of the Milnor fibre. The image M of the counterclock- 
wise generator of vri(T') is called the monodromy operator and fulfills 
M(s)(r) = s(r + 1) for s G if. The sheaf H is determined by the mon- 
odromy representation up to isomorphism. The following well known 
theorem is due to E. Brieskorn [2, 0.6] and others. 

Theorem 1 (Monodromy Theorem). The eigenvalues of the mon- 
odromy are roots of unity and its Jordan blocks have size at most 
(n + 1) X (n + 1) and size at most n x n for eigenvalue 1. 



3. The Gauss-Manin Connection 

Let M = MgMu be the decomposition of M into semisimple part 
M^ and unipotent part M„ and let N = — '°2^" be the nilpotent part 
of M. Note that -27riN G End(^{H(^) where = H (gi^ Q. Let 
He = He be the decomposition of = H ®^ C into generalized 
A-eigenspaces of M and M^ = M|^a. Note that if q = if^^ © ii"^^ 

where ®q C = and H^^ 0q C = ^Xy^i ^c- Then there is an 
inclusion 

defined by ?/'o(/l) = r+^A = exp(X log(t)) with image = im.{ipa)- 
In particular, the operators M and N act on C". The following lemma 
is an immediate consequence of the definition of ipa- 

Lemma 2. 

(1) toip^ = t/^q+i and dtoip^ = ip^_^ o (a + N). 

(2) t : — > C°'~^^ is bijective and dt : — > 0°"^^ is bijective if 

(3) On C", tdt-a = N and exp{-2mtdt) = M^"'"'". 

(4) = ker(t(9t - a)"+^ 

Definition 3. We call G = 0_i<„<o C (i*jr)o the local 

Gauss-Manin connection. 

The local Gauss-Manin connection is a /i-dimensional C{t}[t^^]-vectorspace 
and a regular C{t}[9t]-module. The generalized a-eigenspaces of the 
operator tdt define the decreasing filtration on G by free C{t}-modules 



a<l3<a+l a</3<a+l 
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of rank /x, the V-filtration. In contrast to the ipa and C", the V^" are 
independent of the coordinate t. The define a sphtting 

of the V-fihration and we denote by leady the leading term with respect 
to this sphtting. The ring C{t} is a free module of rank 1 over the ring 

oo oo 

where s = dt acts by integration. This fact is generalized by the 
following lemma ^Sl 1.3.11]. 

Lemma 4. The action of s = d^^ on V^"^ extends to a <C{{s}} -module 
structure and V^~^ is a free <C{{s}} -module of rank fi. 

Since [dt,t] = 1, [t, s] = and hence 

t = s'^ds, dtt = sds- 

We call a free C{{s}}-submodule of V^~^ of rank /i a C{{s}}-lattice 
and call a tc^i-invariant C{{s}}-lattice saturated. A basis e of a C{{s}}- 
lattice defines a matrix A = X]fc>o ^fcS*^ of t by te = eA such that 

t^ A + s^ds 
is the basis representation of t. 

4. The Brieskorn Lattice 

The description of cohomology in terms of holomorphic differen- 
tial forms by the de Rham isomorphism leads to the definition of the 
Brieskorn lattice 

H" = n%'/dfAdnxl 

By E. Brieskorn [21 1.5] and M. Sebastiani ^2], the Brieskorn lattice is 
the stalk at of a locally free ^^-sheaf J^" of rank fi with Jif"\T' = 
and hence H" C {i^J^)o. The regularity of the Gauss-Manin connec- 
tion proved by E. Brieskorn p| 2.2] implies that H" C G. B. Malgrange 
jHJ 4.5] improved this result by the following theorem. 

Theorem 5. H" C V-\ 

By E. Brieskorn |2., 1.5], the Leray residue formula can be used to 
express the action of dt in terms of differential forms by dt [df A uj] = 
[da;]. In particular, sH" C H" and 

H"/sH" = ^]^+Vd/ A QX, - <C{x}/{d{f)). 

Since the V-^'^ is a C{{s}}-module, theorem implies that H" is a free 
C{{s}}-module of rank /i and the action of s can be expressed in terms 
of differential forms by 

s[duj] = [df Aio]. 
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For computational purposes, we may restrict our attention to the com- 
pletion of the Brieskorn lattice. E. Brieskorn [2, 3.4] proved the follow- 
ing theorem. 

Theorem 6. The mx,o- o.'fid xnT,o-adic topologies on H" coincide. 

While the proof of theorem IHl is highly non-trivial, the analogous 
statement for the C{{s}}-structure of the Brieskorn lattice is quite el- 
ementary [ini 1.5.4]. 

Proposition 7. The vnx,o- and vn.c^^s}}-o.dic topologies on H" coincide. 

We call the completion H" of H" the formal Brieskorn lattice. Since 
completion is faithfully flat, H" is a free C|s] -module of rank /i with 
a differential operator t = s'^ds- The equality [d{f)gdx\ = s[d{g)dx\ 
motivates to consider the differential relation d{f) — sd. It is not diffi- 
cult to prove that it defines the formal Brieskorn lattice as a quotient 
of C[s,x] [H 1.5.6]. 

Proposition 8. 

C[s,x] ^ Cls,x}/{d{f) - sd)€ls,x] -cw H". 

Proposition |H1 is the starting point for an algorithmic approach to 
the local Gauss-Manin connection. Let <^ be a local degree ordering 
on C[a;] such that deg(a;) < and deg{d) = — deg(x) > 0. One can 
compute a polynomial standard basis g of the Jacobian ideal {d{f)) and 

a polynomial transformation matrix B = {Vy such that g = d_{f)B. 

By Nakayama's lemma, m = (x^) xi.f(iead{g)) represents a C[s]-basis [m] 

of H". Let <s be the local degree ordering on C|s] and let <= (<s, <^) 
be the block ordering of <s and <^ on C[s,x]. 

Definition 9. 

(1) h= {{g,-sdV)x^)^^^. 

(2) deg(s) = mindeg(m) -|- 2mindeg(x) < 0. 

(3) = {Nk)k>o with Nk = K deg{s) - 2mindeg(x). 

(4) V = {Vk)k>o with VK = {pe Cls,a;]| deg(p) < Nk} + (s)^ C 
€ls,xl 

Since H" is a free C|s]-module, h is a. standard basis of the C[s]- 
module {d{f) — s9)C[s,a;]. The following lemma is technical but not 
very deep and can be generalized to formal differential deformations 
[ini 2.2.10]. 

Lemma 10. V = {Vk)k>o is a basis of the {s,x)-adic topology of 
€ls,x] with TTHiVK) = {s)^H". //s"lead(/ij-^) G Vk then s'^h^^p e Vk- 

Lemma ITUl leads to a normal form algorithm for the Brieskorn lattice 
[T^ 2.2.12]. It computes a normal form with respect to h and hence 
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the [m]-basis representation in H" . The normal form computation up 
to a given degree can be continued up to any higher degree without 
additional computational effort. The normal form algorithm for the 
Brieskorn lattice is a special case of a modification of Buchberger's 
normal form algorithm jSj for power series rings where termination is 
replaced by adic convergence [ini 2.1.19]. 

5. Mixed Hodge Structure 
By lemma 121 there is a C-isomorphism 

-l<a<0 -l<a<0 

defined hj ip = ©-i<q<o "^a and the monodromy M on Hp, corresponds 
to exp{-2mtdt) on 0_i<„<oC'"- 

The Hodge filtration F = {Fk)k&'z on V^^^ defined by J. Scherk and 
J.H.M. Steenbrink is the increasing filtration by the free C{{s}}- 
modules 

Fk = F"-'^ = {s^'^H") n V>-^ 

of rank fi. Via the splitting C" = giyV^'^, the Hodge filtration induces 
an increasing Hodge filtration FC"' by C-vectorspaces on and, via 
ip, on Hp. The nilpotent operator — 27riN G EndQ(ifQ) defines an 
increasing weight filtration W = {Wk)ke'z centered at n resp. n + 1 on 
Ht^ resp. H^. 



Theorem 11. The weight filtration W on Hq and the Hodge filtra- 
tion F on Hp define a mixed Hodge structure on the cohomology H 
of the Milnor fibre and the operator N is a morphism of mixed Hodge 
structures of type (—1,-1). 

The mixed Hodge structure on the cohomology of the Milnor fibre 
was discovered by J.H.M. Steenbrink and described in terms of the 
Brieskorn lattice by A.N. Varchenko jTTI . 

The nilpotent operator on defines an increasing weight filtra- 
tion W = (Wk) k^w. centered at n on C". By definition N commutes 
with ipa and hence 



W[-1]C^, aeZ. 



The weight filtration W = 0_i<„<o C{{s}}VrC° on V>-^ by free 
C{{s}}-modules induces WC^ via the splitting = gr^l^^"^. 

The spectral pairs are those pairs (a, /) G (Q x Z with positive mul- 
tiplicity 

rfr = dimcgrrg4gr^^>-\ 
Via the isomorphism ip, they correspond to the Hodge numbers 

hP/~P = dime gi^pgr^H^ 
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by d'^+P = for -1< a < and < = and inherit 

the symmetry properties 

ja j2n—l—l—a ja ja—n+l ja jn—l—a 

from the mixed Hodge structure. The spectral numbers are those num- 
bers a G Q with positive muhiphcity 

d'' = dimegr'^gr^V>-' = J2d? 

lew. 

and have the symmetry property d" = 

6. M. Saito's Basis 

By P. Dehgne j^J 1.2.8], a morphism of mixed Hodge structures is 
strict for the Hodge filtration. In particular, by theorem ITU N is strict 
for the Hodge filtration on and on gr^V^-^"^. Hence, there is a direct 
sum decomposition F^C" = ^.^^C^^i such that N(C"''=) c C"'*^+i, 

and sC"''^ C By definition of the Hodge filtration, 

leadv(i7") = J2I1 'C{{s}}C-'>^ = €{{s}}G- 

where = C"'°. Let <qxz= (>Q) >z) be the block ordering of 
>Q and >z on the index set Q x Z. Then the Hodge filtration de- 
fines a refinement of the V-filtration on by free C{{s}}-modules 
V^^^ = FkC^ © 1/>" of rank /i and the C''' define a splitting of this 
refined filtration compatible with s. We call the refinement the Hodge 
refinement and the splitting a Hodge splitting. The following lemma 
follows essentially from the fact that C{{s}} is a discrete valuation ring 
[131 1.10.5,1.10.10]. 

Lemma 12. Let H be a C{{s}} -lattice and C"''^ a splitting of a refine- 
ment of the V-filtration compatible with s. Then a minimal standard 
basis of H is a C{{s}}-&aszs and there is a reduced minimal standard 
basis of H . 

In particular, there is a reduced minimal standard basis of H" for 
a Hodge splitting. The following proposition follows essentially from 
lemmaElS [13, 1.10.12]. 

Proposition 13. Let h be a reduced minimal standard basis of H" for 
a Hodge splitting. Then the h-matrix A oft has degree 1. In particular, 

{H", t) ^ {€{{s}r, Ao + A^s + s'ds) 

is an isomorphism. Moreover, Ai is semisimple with eigenvalues the 
spectral numbers of f added by 1 and giy{AQ) can be identified with N. 
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Note that the matrices Aq and Ai in proposition determine the 
differential structure of the Brieskorn lattice. M. Saito [lOj first con- 
structed a C{{s}}-basis of H" as in proposition IT^ without calhng it a 
reduced minimal standard basis. 

7. The Algorithm 

We describe an algorithm to compute Aq and Ai as in proposition 
ITHl This algorithm can be simplified to compute the complex 

monodromy, the spectral numbers, or the spectral pairs only [T!^ . 

The normal form algorithm for the Brieskorn lattice in section |3] 
computes the [m]-matrix A = Ylk>o ^kS^ of t defined by t[rn\ = [frn\ = 
[rn\A up to any degree. We identify the columns of a matrix H with 
the generators of a submodule {H) C Cls^ and denote by E the unit 
matrix. Then (E) is the [m]-basis representation of H". Hence, the 
following two statements hold for h = [m] with n = and H = E. 

{Hh) One can compute n > and a /i x /x-matrix H with coeffi- 
cients in C[s] of degree at most k such that (H) is the /i-basis 
representation of H" and s'^{E) C (H). 
(Ah) One can compute the /i-matrix A of t up to any degree. 
Step by step, we improve the C[s]-basis h and show that (Hh) and 
(Ah) hold. After the last step, Aq and Ai as in proposition can 
be computed by a basis transformation of A to a reduced minimal 
standard basis of {H) up to a certain degree bound. 

We call the canonical projection jet^ : C[s] — > 0^=0^5-' the k- 
jet. Let the monomial ordering on Cls]'^ = C|s] (S>c C'^ be the block 
ordering <= (<s, >^) of the local degree ordering <s on C|,s] and the 
inverse ordering >^ on the indices of the basis elements of C^. 

7.1. The Saturation of H". In this step, we show that {Hh) and {Ah) 
hold for a C|s]-basis hof a saturated C[s]-lattice. 
The increasing sequence of C|s]-lattices defined by 

H'o = H", H'l^^ = sH'l + tH'l C H" 

is stationary since H" is noetherian. Hence, the saturation H'^ = 
lJj^.>Q H'l^ of H" is a saturated C[s]-lattice. The [m]-basis representation 
{Hk) of H'f! can be computed by 

Ho = g_i = E, Qk= (jetfc(A) + s^ds)Qk-u H^+i = {sHk\Qk). 

We successively compute the and check in each step if {Qk) C {Hk) 
by a standard basis and normal form computation. If {Qk) C {H^) 
then we stop the computation and set k = and H^o = H^. Then 
{Hoo) is the [m]-basis representation of H'^. We replace by a 
minimal standard basis of (Hoo). Then h = s~'^hHoo is a C[s]-basis of 
a saturated C[s]-lattice. By a normal form computation with respect to 
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Hoc up to degree k, we compute the /i-basis representation (H^s'^E) = 
iiet^iH-h'^E)) oiH". Since {H^) C (E), s-{E) C {H-^s^E). By a 
normal form computation with respect to H^o up to degree k, + k, one 
can compute the fc-jet 

ietj,{H-\A-KsE+s^ds)H^) = jet^{H-' {iet,^,,{A-KsE)+s^ds)H^) 
of the /i-matrix of t for any A; > 0. 

7.2. The V-Filtration. In this step, we show that (if^) and (A^) hold 
for a <Q-increasingly ordered C|s]-basis v of a V°' compatible with the 
direct sum decomposition V"' /sV"' = ®a<(3<a+i 

Since his a C[s]-basis of a saturated C{{s}}-lattice, Aq = and, by 
theorem the eigenvalues of Ai are rational. In order to compute the 
eigenvalues of Ai, we transform Ai to Hessenberg form and factorize 
the characteristic polynomials of its blocks. Then we compute a con- 
stant C[s]-basis transformation such that Ai = diag(Q;i, . . . , a^) + N 
with ai < ■ ■ ■ < where diag(ai, . . . , a^) denotes the diagonal matrix 
with entries ai, . . . ,a^. If — ai < 1 then v = hisa <Q-increasingly 
ordered C[s]-basis w of a V"' compatible with the direct sum decom- 
position V'^/sV"' = ®a<f3<a+i ■ If ttfi — tti > 1 then we proceed as 
follows. Let 



such that Af) = 0, A}'^ = 0, A^'^ = 0, and the eigenvalues of A}'^ are 
the eigenvalues a of Ai with a < ai + 1. Then the C|s] [s^^]-basis 
transformation 

decreases a^ — ai and the degree up to which A is computed by 1 and 
increases k by 1. After at most n such transformations, — ai < 1. 

7.3. The Canonical V-Splitting. In this step, we show that {He) 
and (Ac) hold for a <Q-increasingly ordered C-basis c of a direct sum 

a</3<a+i compatible with the direct sum. 

Let c be the image of [y] under the splitting V"/ sV" = ®a<f3<a+i 
By Nakayama's lemma, c is a C-basis of @a<i3<a+i compatible with 
the direct sum. The eigenvalues of the commutator [-, Ai] G Endc(C^^) 
are the differences of the eigenvalues of Ai. Since — ai < 1, 
- k e GL^2(C) for k >l. Let f/ = Ej°lo^i*-' be the C[s]- 
basis transformation defined by c = yJJ . Then Uq = E and UAis = 
{A + s'^ds)U or equivalently 

fc-i 

Uk = {[■,Ai]- ky^J2^k-j+iUj 

j=0 
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for A; > 1 and hence one can compute U up to any degree. Since Uq = E 
and K > 0, jet^(?7) is a minimal standard basis of (E). By a normal 
form computation with respect to U up to degree k, we compute the 
c-basis representation {U^^H) = (jet^(jet^(f/)^^if)) of H" and Ais is 
the c-matrix of t. 

7.4. A Hodge Splitting. 

In this step, we show that (Hf) and (Af) hold for a <Qxz-decreasingly 

ordered C-basis / of a direct sum ®a<f3<a+i ©fcez ^'^''^ compatible 
with the direct sum and that one can compute Aq and Ai as in propo- 
sition uni 

We compute a standard basis of H up to degree k in order to com- 
pute the c-basis representation of the Hodge filtration F. The nilpotent 
part of Ai is the c-basis representation of N. By computing images and 
quotients of C-vectorspaces, we compute the c-basis representation of a 
Hodge splitting FkC^ = 0j<^ C^'^ . Then we compute a constant C[s]- 
basis transformation / = cU such that / is a <Qxz-decreasingly or- 
dered C-basis of the direct sum ®a<f3<a+i ©itez compatible with 
the direct sum. 

We replace if by a reduced minimal standard basis of (H) up to 
degree k + 1. By a normal form computation with respect to H up to 
degree n + 1, we compute the 1-jet 

iet,{H-\A + s'ds)H) =jeti(jet,+i(iJ)-i(jet,+i(A) + s29,)jet,+i(i7)) 
of the cif-matrix A of t in order to compute Aq and Ai as in proposition 

cni 

8. An Example 

The algorithm in section [3 is implemented in the computer alge- 
bra system Singular [0] in the procedure tmatrix in the library 
gaussman. lib ^21- In an example SINGULAR session, we compute 
the differential structure of the Brieskorn lattice of the singularity of 
type T2,5,5 defined by the polynomial / = x'^y'^ + + y^. 

First, we load the Singular library gaussman. lib: 

> LIB "gaussman. lib" ; 

Then, we define the local ring R = Q,[x,y]{x,y) with the local degree 
ordering ds as monomial ordering and the polynomial / = x^y^ + x^ + 
G R: 

> ring R=0, (x,y) ,ds; 

> poly f=x2y2+x5+y5; 

Finally, we compute Aq and Ai as in proposition UHl 

> list A=tmatrix(f ) ; 
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The result is the hst A=A[1] ,A[2] such that A[i + 1] = Ai and 
/O ••• o\ 



An 



••• 

V ••• 0/ 



/1^^_9__9_111113133 
Ai - diag(^-, — , — , — , — , 1, — , — , — , — , - 



By proposition Uni {H",t) ^ (C{{s}}^,Ao + sAi + s^d.,) and the 
spectral pairs are (-|,2), (-3^, l)^ (0, 1), 

(io). 
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